arXiv:1508.06794vl [astro-ph.HE] 27 Aug 2015 


PROCEEDINGS 

0F SCIENCE 



The search for short-term flares in extended VHE 
Crab Nebula observations with the Whipple 10 m 
telescope 


A. O’Faolain de Bhroithe* for the VERITAS collaboration 1 

DESY, Platanenallee 6, 15738 Zeuthen, Germany 
E-mail: anna . of aolain . de . bhroithe@desy. de 

In 1989, the Whipple 10 m telescope achieved the first indisputable detection of a TeV gamma-ray 
source, the Crab Nebula. Until its decommissioning in 2011, the Whipple telescope took regular 
measurements of the nebula. With the recent discovery of GeV gamma-ray flaring activity in 
the Crab Nebula, it is an opportune time to return to the Whipple telescope data set and search 
its extensive archive for evidence of TeV flares. A data set on the Crab Nebula spanning ten 
years, 2000 - 2010, is compiled and searched for day-scale flaring activity using a Bayesian-block 
binning algorithm. No evidence for significant flaring activity is found. Monte Carlo simulations 
show that low levels of flux increase on short timescales are difficult to detect. Assuming a 
flare duration of seven days, 99% confidence level upper limits are calculated for the possible 
frequency of five-fold, two-fold and 1.5-fold flares in the data set. An upper limit of 0.02 flares 
per year is found for the five-fold flare, and a limit of 0.27 flares per year is placed on the two-fold 
flare. The detection of the 1.5-fold flare is consistent with the false-positive rate of the method, 
and so cannot be excluded. 


The 34th International Cosmic Ray Conference, 
30 July- 6 August, 2015 
The Hague, The Netherlands 


* Speaker. 

Neritas. sao.arizona.edu 


© Copyright owned by the author(s) under the terms of the Creative Commons Attribution-NonCommercial-ShareAlike Licence. 


http://pos.sissa.it/ 








Crab Nebula observations with the Whipple 10 m telescope 


A. O’Faolain de Bhroithe 


1. Introduction 

The supernova explosion that created the Crab Nebula and its central power source, a pulsar, 
was first observed in 1054 AD. The Crab Nebula was the first source detected at TeV energies 
[1] and its discovery laid the foundation for the field of very-high-energy (VHE; E > 100 GeV) 
astrophysics. It has been considered a standard VHE reference source, similar to its role in X-ray 
and high-energy gamma-ray astronomy, e.g., [2, 3]. 

In the late 1960s, radio and optical pulsations were detected from what is now known to be 
the Crab pulsar [4, 5] with a spin-down rate of 36 ns per day [6]. A magnetohydrodynamic model 
was developed using the pulsar as the energy source of the nebula [7-9]. This model was capable 
of accounting for most of the observed features (see Figure 1 for a schematic illustration) and 
predicted constant emission from the nebula. 

In 2011, [10] reported a decline of ~ 7% in the Crab Nebula flux in the 15-50keV band 
from 2008 to 2010 observed with the Fermi-GBM and independently confirmed with the Swift- 
BAT, the RXTE-PCA and the INTEGRAL-IBIS. The same publication reported a similar decline 
in the 3-15 keV band (RXTE-PCA), and in the 50- 100 keV band ( Fermi-GBM , Swift-BAT and 
INTEGRAL-IBIS). Contemporaneous with the publication of the X-ray variability, the discovery of 
short-term flares on the timescale of days with the Fermi -LAT [11] and AGILE [12] was announced. 
The Fermi-LAT team reported > 100 MeV flares detected in 2009 February and 2010 September, 
and the AGILE team reported flares in the range 0.1 - 10 GeV in 2007 October and 2010 September, 
the latter being the same flare as that observed with the Fenni- LAT. In April of 2011, the Fermi- 
LAT detected the strongest gamma-ray flare seen to date [13]. This flare lasted approximately nine 
days and reached a peak flux corresponding to a 30-fold increase in flux above 100 MeV. 

The discovery of large flux increases on timescales of days and shorter was completely un¬ 
expected in the framework of the existing theoretical model. New models have been proposed to 
explain the newfound behavior, often suggesting regions of the nebula corresponding to the optical 
and X-ray features as the acceleration sites. In general, the models suggest one of three scenarios: 
instabilities in the termination shock, i.e., the inner knot feature e.g., [14], particle acceleration 
due to magnetic reconnection events in the nebula e.g., [15], or plasma kink instabilities e.g., [16] 
which may provide a framework for variability originating in the anvil region. Of these, only the 
model of [15] provides predictions for the VHE domain. 

Both MAGIC and VERITAS reported no change in flux or spectral shape of the VHE Crab 
Nebula during the GeV flare in 2010 September [17, 18]. However, ARGO-YBJ reported a flux 
about 3-4 times higher than the average emission at a median energy of about 1 TeV with a sta¬ 
tistical significance of ~ 4 a, based on a preliminary analysis of their data [19]. A more recent 
publication sets the posttrials significance of this increase to be of the order of ~ 2a [20]. This 
publication also reports the ARGO-YBJ long-term light curves of the Crab Nebula to be consistent 
with a uniform flux for median energies of 0.76, 1.8, and5.1 TeV respectively. Even though the rate 
variations in the ARGO-YBJ light curves appear to be consistent with statistical fluctuations, a cor¬ 
relation study of the variation in the event rate measured by ARGO-YBJ (median energy 0.76 TeV) 
and the flux measured by the Fermi-LAI (above 100 MeV) shows a Pearson correlation coeffi¬ 
cient of 0.56 ±0.22, which could hint at similar long-term behavior of the gamma-ray emission at 
~100MeV and ~1 TeV. 
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Figure 1: Schematic illustration of the characteristic features of the Crab nebula seen at optical and X-ray 
energies. 


During another GeV gamma-ray flare from the Crab Nebula in 2013 March, in which a 20- 
fold increase in flux above 100 MeV for the synchrotron component of the nebula was detected 
[21], simultaneous VHE observations were taken by VERITAS [22] and H.E.S.S. [23]. Neither 
instrument detected evidence of any VHE flux enhancement. Since the discovery of the flares, 
many other simultaneous multiwavelength observations across the electromagnetic spectrum have 
been obtained in order to search for correlated emission at different energies, but no such correlation 
has been found e.g., [13, 24-27], 

Coverage of the GeV flares with the Whipple 10 m telescope is unfortunately sparse, with only 
two days of overlap at the very end of the 2007 October flare (showing nothing unusual) and no 
overlap with the 2009 February flare. However, it is reasonable to assume that there have been 
flares from the nebula prior to the current gamma-ray satellite era. The Whipple telescope has 
a large archive of VHE data on the Crab Nebula dating back to the initial detection of the Crab 
Nebula with the telescope in 1989. In this work, a ten-year data set collected with the Whipple 
10 m telescope on the Crab Nebula from 2000 to 2010 is searched for flaring activity. Throughout 
this period, the configuration of the Whipple telescope remained relatively consistent from year to 
year, making it particularly suitable for variability studies. 

The 2000-2010 Whipple telescope observations are described in § 2. A Bayesian-block bin¬ 
ning algorithm is used to search for variability in the data set, and is presented along with the results 
in § 3. A discussion of these results follows in § 4. 

2. Observations and event reconstruction 

The Whipple 10 m telescope [28] was the first large, purpose-built telescope for atmospheric - 
Cherenkov astronomy. Located at the Fred Lawrence Whipple Observatory in southern Arizona 
at an altitude of ~2.3km above sea level, it was completed in 1968 and operated as an imaging 
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atmospheric-Cherenkov telescope (IACT), i.e., with an imaging camera, from 1982 to 2011, after 
which it was decommissioned. It was sensitive to photons with energies up to 20TeV, with a 
nominal low-energy threshold (which remained constant over the ten-year period of interest) of 
~ 400 GeV for a Crab Nebula-like spectrum. 

The telescope was of Davies-Cotton design, comprising a 10 m diameter reflector with a ra¬ 
dius of curvature of 7.3 m and the camera at the focal plane. The reflector was composed of 248 
tessellated hexagonal mirrors with a radius of curvature twice that of the overall dish. The camera 
was upgraded many times over the lifetime of the telescope. The two camera configurations used 
in this work each had a trigger field of view of ~ 2.6°. The different cameras and their properties 
are summarized in [28]. 

The telescope system was shut down between the end of June and the start of September 
each year to protect the electronic equipment from lightning damage due to the monsoon season 
in Arizona. As a result, an observing “season” is defined as beginning in September and lasting 
approximately until June. The Crab Nebula is visible at zenith angles < 35° from the Whipple 
Observatory site between August and April each year, allowing good quality observations of the 
source for up to seven months each observing season. The telescope did not operate for seven 
nights around the full moon period, resulting in a seven-day gap in each lunar cycle (~28 days). 

A data set was compiled using Crab Nebula observations taken with the Whipple 10 m tele¬ 
scope during the years 2000-2010. Data were taken in 28-minute observations in paired mode , in 
which observations of the source were paired with observations of a region of the sky at the same 
azimuth and elevation angles with no known gamma-ray emitter for the purposes of background 
estimation. Data were required to be accumulated under good weather conditions and at zenith 
angles < 35°, where the effective area is expected to be approximately constant. A total of ~ 150 
hours of quality-selected live time satisfied these criteria, including 328 observation pairs taken on 
243 separate nights. 

The data were analyzed using the standard Supercuts procedure described in Appendix B of 
[29]. This analysis is a two-step process. In the first step, the images arc processed to remove hard¬ 
ware and night-sky-background effects. This is achieved by flatfielding the images using relative- 
gain information from nightly calibration runs taken with a nitrogen arc lamp, screened with a 
diffuser to ensure uniform illumination of the camera. Software padding is used to compensate for 
any difference in sky brightness between the source and background regions, by injecting software 
noise into the events for the darker region [30]. In the second step, the moment-fitted image pa¬ 
rameters are calculated [31] and the candidate gamma-ray events selected based on the principal 
components and derived parameters of the images. 

A light curve of the Crab Nebula, binned by individual 28-minute observation and expressed 
in units of gamma rays per minute, is created for each observing season. The light curves arc 
separated by season as no interseason calibration has been applied. Each light curve is fitted with a 
constant rate to test for temporal stability. The rate for each season is found to be consistent with a 
constant fit within 3 a. 

3. Bayesian-block analysis 

The data set is searched for variability using the Bayesian-block binning algorithm [32]. This 
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is a nonparametric modeling technique that is applied to the general problem of detecting and 
characterizing local variability in sequential data. The method is used to obtain a representation of 
the light curve in terms of a segmentation of the time axis into blocks — subintervals of generally 
unequal size containing consecutive data that satisfies some well-defined criterion. A good example 
of such a criterion is the sum over all blocks of a goodness-of-fit measure of a given model of the 
data in each block. The optimal segmentation is that which maximizes this goodness-of-fit. In 
this analysis, the piecewise-constant model is used to describe the data in each block. This is the 
simplest model and allows for an exact treatment of the likelihood used to calculate the overall 
goodness-of-fit. 

The form of the likelihood differs for different types of data. The type of data used in this 
analysis — source rates measured at different times — is categorized as point measurements ac¬ 
cording to [32, §1.5], and the corresponding form of the likelihood is used. A false-positive rate 
of po = 0.01 is selected for this analysis as a compromise between a low false-positive rate and 
retaining the sensitivity of the test to fluctuations in the signal. 

Figure 2 shows the ten-season light curve of the Crab Nebula, binned by individual 28-minute 
observation. The result of the Bayesian-block binning on the data set is also shown. The light curve 
is best represented by a single time bin spanning the entire data set, which implies that the data is 
consistent with a constant source emission model. 

4. Discussion 

No evidence for variability is found in ten years of Whipple telescope Crab Nebula observa¬ 
tions. However, it is possible to obtain a limit on the level of flaring activity that could be detected. 
In order to do this, a toy Monte Carlo simulation was developed to simulate a single flare of known 
length and emission within an otherwise standard data set. Simulations of observations were cre¬ 
ated using separate Poisson distributions with constant mean for the source and background rates. 
The means of these distributions were set to those observed in the real data. No observations were 
simulated for seven consecutive nights out of each 28, mimicking the bright moonlight period. For 
the remaining nights, the number of Crab Nebula observations on each night was randomly chosen 
from an empirical distribution of the source sampling obtained by counting the number of Crab 
Nebula observations per night in the real Whipple telescope data (numbers range from zero to four 
observations of the source per night). A flare was then injected beginning on a random day, with 
the possibility that flares may overlap the bright moonlight period when they cannot be observed. 
A medium-duration flare of seven days was used for this study. It was found that on average, 
~30% of the injected flares go unobserved due to a combination of the bright period and the data 
sampling. 

To quantify the performance of the Bayesian-block algorithm, the simulation was run 10000 
times each for flare emission levels corresponding to five times, two times, and 1.5 times the av¬ 
erage Crab Nebula flux. The five-fold flare was detected in all simulations with at least one obser¬ 
vation during the flare period. Using the method of [33], a 99% confidence level (CL) upper limit 
was calculated on the frequency of flares that may be present in the data set. An upper limit of 0.02 
five-fold flares per year was found, assuming at least a single observation of the flare. The two-fold 
flare was detected in 37% of the simulations with at least one observation of the flare, leading to a 
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Figure 2: Ten-year light curve of the Crab Nebula covering the seasons 2000-2010, binned by individual 
28-minute observation. The optimal temporal binning determined by the Bayesian-block analysis is a single 
bin spanning the entire data set, and is indicated by the solid gray line. 


99% CL upper limit, of 0.27 flares per year of this magnitude. The 1.5-fold flare was detected in 
1% of cases, consistent with the false-positive rate of this analysis, and so cannot be excluded. 

The only model in the literature to date with predictions for the TeV behavior of the source 
in the context of the flares is provided by [15]. These authors suggest that the variable gamma- 
ray emission originates behind the shock where the pulsar wind interacts with the nebula as a 
result of magnetic reconnection. The model predicts TeV variability on the same timescales as 
those observed from the synchrotron nebula, as it is the same population of synchrotron-radiating 
electrons that upscatter soft photons into the VHE regime. Such variability is expected to be of the 
order of ~10% above 1 TeV, and more substantial above 10 TeV. As it is difficult to detect flares 
of the order of 50% flux increase, a limit cannot be placed on the frequency of low level of flaring 
activity that may be present in the source according to this particular model. 

The nondetection of any flaring behavior in the Whipple 10 m telescope data set is consistent 
with the VERITAS [22] and H.E.S.S. [23] observations of the 2013 March GeV 11 are, which found 
no enhancement of the VHE flux on the level of a few percent. The detection of TeV variability 
in the Crab Nebula would be a very important result for determining the fl are mechanism. It is 
possible that future instruments, especially those with a larger collection area for TeV gamma rays 
such as the Cherenkov Telescope Array, will be able to detect small variations in the VHE flux 
from the Crab Nebula on short timescales [34], or place even stronger constraints on the level of 
VHE variability of the nebula during the high-energy synchrotron flares. 
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